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Marginally trapped surfaces (MTSs) are commonly used in numerical relativity to locate black 
holes. For dynamical black holes, it is not known generally if this procedure is sufficiently reliable. 
Even for Schwarzschild black holes, Wald and Iyer constructed foliations which come arbitrarily close 
to the singularity but do not contain any MTSs. In this paper, we review the Wald-Iyer construction, 
discuss some implications for numerical relativity, and generalize to the well known Vaidya space- 
time describing spherically symmetric collapse of null dust. In the Vaidya spacetime, we numerically 
locate non-spherically symmetric trapped surfaces which extend outside the standard spherically 
symmetric trapping horizon. This shows that MTSs are common in this spacetime and that the event 
horizon is the most likely candidate for the boundary of the trapped region. 
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Introduction: In stationary black hole spacetimes, there 
is a strong correspondence between marginally trapped 
surfaces (MTSs) and event horizons (EHs) because cross 
sections of stationary EHs are MTSs. MTSs also feature 
prominently in the frameworks of isolated, dynamical 
and trapping horizons which have shed considerable 
light on the properties of classical and quantum black 
holes even in the dynamical regime. See e.g. flUi] 
for reviews. Numerical simulations routinely look for 
MTSs to locate black holes on a Cauchy surface. This is 
because, while MTSs can be located on a Cauchy surface 
in real time, the EH can only be located a posteriori after 
the simulation has been successfully completed. MTSs 
can be useful for extracting physical information about 
a black hole in a numerical simulation 0]. 

However, in dynamical situations, the correspon- 
dence between MTSs and EHs is lost (beyond the fact 
that MTSs are enclosed by the EH); the event horizon is, 
in general, an expanding null surface, while outgoing 
light rays from a MTS have, by definition, zero expan- 
sion. An explicit example was constructed by Wald and 
Iyer [5] where they showed that even in the Schwarz- 
schild spacetime, there exist perfectly regular Cauchy 
surfaces which come arbitrarily close to the singularity 
and foliate the spacetime, but which nevertheless do not 
contain any MTSs. While this has not been an issue in 
most numerical simulations to date, it raises the ques- 
tion of whether MTSs can be found generally in numer- 
ical simulations of black hole spacetimes. 

However, there are other results which indicate that 
MTSs should be common in black hole spacetimes. For 
example, it was suggested by Eardley 0|, that an MTS 
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can be locally perturbed in an arbitrary spacelike direc- 
tion to yield a 1-parameter family of MTSs. A precise 
formulation of this idea and its proof follows from re- 
cent results of Andersson et al. |7]. From this result, it 
seems plausible that for a "generic" Cauchy surface S 
passing through the black hole region, one should be 
able to perturb a nearby MTS to lie on Z,. Thus, Cauchy 
surfaces passing through the black hole should "gener- 
ically" contain a MTS. A related issue is the question 
of the boundary of the trapped region of a black hole 
spacetime; clearly, a MTS cannot be perturbed to lie out- 
side this boundary. Eardley suggests that the boundary 
should be the event horizon, while arguments by Hay- 
ward suggest that the boundary should be a trapping 
horizon |8]. For stationary black holes the two notions 
coincide, but not for dynamical black holes. 

In this paper, we study these issues from a numerical 
relativity perspective. In particular, we study the spher- 
ically symmetric Schwarzschild and Vaidya spacetimes 
using foliations not adapted to the spherical symmetry. 
Very little is known so far about trapped or marginally 
trapped surfaces on such slices, either analytically or nu- 
merically. We find that for Vaidya, MTSs are indeed 
easy to locate, and we do not encounter the problem 
suggested by Wald and Iyer. This suggests that Cauchy 
surfaces of the type suggested by Wald and Iyer (which 
presumably exist also in the Vaidya spacetime) are ex- 
ceptional. In Vaidya, we find that the non-symmetric 
marginal surfaces lie partially outside the standard r = 
2M surface (which in this case is a trapping horizon, 
and lies inside the event horizon), thus indicating that 
the event horizon is a better candidate for the bound- 
ary of the trapped region as suggested by Eardley. In 
the remainder of this paper, we review basic concepts 
regarding trapped surfaces and horizons, outline Wald 
and Iyer 's construction with an explicit example and re- 
mark on its implications for numerical relativity, and fi- 
nally discuss trapped surfaces in the Vaidya spacetime. 
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Trapped surfaces and horizons: Let l" and n" be the 

future directed null-normals of a closed 2-surface S. Let 
qah be the 2-metric on S induced by the spacetime metric. 
The expansion of F is 0(^) = q°'^Vah with a similar 
definition for 0(„). The surface S is said to be trapped 
if both expansions are negative: < and &(n) < 

0. For a marginally trapped surface (MTS), one or both of 
these inequalities are replaced by an equality instead. 
Weakly trapped surfaces have 0(„) < 0, < 0. All 
these definitions are invariant under arbitrary positive 
rescalings of £" and n". 

The definition of a marginally outer-trapped surface 
(MOTS) requires a choice of an "outgoing" direction 
with respect to future null infinity or spatial infinity. 
This choice of outgoing direction breaks the symmetry 
between the two null normals i" and n". A MOTS is thus 
a MTS with 0(f-| — 0, where £ is the outgoing direction. 

The trapped region is the region where trapped surfaces 
exist. This is defined either in the full spacetime or on 
a Cauchy surface E. A point is in the trapped region 
if there is a trapped surface which contains that point. 
Similarly, a point is in the trapped region of 'L iff there is 
a trapped surface on E that contains this point. The ap- 
parent horizon (AH) on Z is the boundary of the trapped 
region of E. As such, its definition is so complicated that 
it is numerically not feasible to look for it directly. How- 
ever, an AH is also a MTS 1 9], and these can be efficiently 
detected. 

Finally, a marginal surface (MS) [8] is a surface where 
one of the null normal's expansion vanishes, i.e., 0(f) = 
where £ can be any of the two null directions, with 
no restriction on 0(,i). Unlike the definition of a MOTS, 
marginal surfaces do not require globally defined out- 
going/ingoing directions. It is called a future marginal 
surface if 0(„) < and a past marginal surface if 0(„) > 0. 
A future marginal surface is the same as a MTS and usu- 
ally arises in numerical simulations as the cross-section 
of a dynamical horizon (DH) 1 10] or an isolated horizon (IH) 
1 4], or more generally, a trapping horizon |8]. 

The Wald-Iyer construction: Wald and Iyer |5] con- 
struct a foliation of the extended Schwarzschild space- 
time in which the spacelike hypersurfaces come arbi- 
trarily close to the singularity, but nevertheless do not 
contain any trapped surfaces. It should be noted that 
these foliations, while special, are not pathological in 
any sense and they can be readily constructed in a nu- 
merical code. Wald and Iyer prove that, if the intersec- 
tion of the slice with the trapped region lies in the past 
of, roughly speaking, "a single event on the future sin- 
gularity", no slice of such a foliation contains a trapped 
surface. This construction relies on the existence of an- 
gular horizons in the black hole region, just as in a cos- 
mological spacetime near the initial singularity. 

An explicit example of one such Cauchy surface is 
easy to construct. Consider the extended Schwarzschild 
spacetime in Kruskal coordinates {T,X,6,<p) (see eq. 
(6.4.29) of |11]). The hypersurface T — k cos can be eas- 



ily shown to satisfy the Wald-Iyer condition for k < 1/2. 
Thus, even though this surface enters the black hole re- 
gion, it does not contain any trapped surfaces. Such a 
slice is depicted schematically in figure ^ It intersects 
the black hole horizon for T > 0, and the white hole 
horizon for T < 0. 

Even though it is a fact that the above Cauchy surface 
does not contain a MOTS, standard apparent horizon 
trackers employed in numerical simulations will hap- 
pily find an "apparent horizon" on this slice. This appar- 
ent contradiction is an issue of terminology. What the 
apparent horizon tracker will locate is the intersection 
of the Cauchy surface with the surface T = X, which is 
the bold line in figure The intersection is the 2-sphere 
given by T = X = fc cos 6. In numerical relativity, one 
typically chooses that part of X+ which belongs to one 
specific asymptotically flat end of the spacetime. Thus, 
the "outgoing" null normal £" and the ingoing null nor- 
mal n" are the ones shown in figure^ With this choice of 
£", the surface given above satisfies 0((!) ~ 0. However, 
this "apparent horizon" is not a MTS because 0(„) < 
on the black hole portion and 0(„) > on the white hole 
portion. 

What is often loosely called "apparent horizon" in 
numerical relativity, or almost as loosely "marginally 
outer-trapped surface", is really only a marginal surface, 
or a future marginal surface if the condition O^,,) < is 
checked (which it is often not). Determining the glob- 
ally outgoing direction is usually either unpractical or 
impossible in numerical simulations. If done, it requires 
some additional knowledge of the simulated spacetime 
that the code itself generally does not have. The easiest 
way to avoid such situations in numerical relativity is to 
explicitly make sure that the apparent horizon is future 
trapped by verifying that &(n) is negative everjrwhere, 
both in the initial data and during evolution. Regarding 
black hole initial data, the construction presented in 1 1^ 
will ensure that the marginal surfaces are future trapped 
and will therefore avoid any of the Wald-Iyer slices if 
the lapse function is kept nonnegative everjrwhere. 

Trapped surfaces in the Vaidya spacetime: Generaliz- 
ing to dynamical situations, consider the Vaidya space- 
time which describes spherically symmetric collapse of 
null dust (radiation) lisll . This is an astrophysically un- 
realistic toy model, but it does serve as a very useful test- 
ing ground. It has been extensively used, for example, 
to study the formation of naked singularities. In ingoing 
Eddington-Finkelstein coordinates {v, r, 6, cp), the metric 
is 

ds^ = - (^l-^^^^^dv^ + ldvdr + r^dCi^, (1) 

where the mass function M{v) can be specified as a 
function of the null coordinate v. For constant M{v), 
this is just the standard Schwarzschild metric in ingo- 
ing Eddington-Finkelstein coordinates. The stress en- 



Figure 1: Penrose diagram of the extended Schwarzschild 
spacetime, from a numerical relativity point of view. The bold 
line is the null world tube of marginal surfaces. The hatched 
area is the region of spacetime that is invisible to the "observer 
at infinity", which is located at the X+ to the right. Every point 
in this figure is a sphere as in conventional Penrose diagrams, 
except for the curves labeled NP and SP. If we represent the 
Cauchy surface as f{T,X, 9, ip) = 0, then NP and SP are re- 
spectively the projections of the north pole {6 = 0) and south 
pole (0 = tt); the intermediate angles lie in between. NP enters 
the trapped region, but SP does not. 

ergy tensor is determined by the derivative of M{v): 

where M = dM/dv > 0. We shall use a time coordinate 
defined as f = c — r and we shall take the mass fimction 
to be non-zero only for c > 0. Thus, the spacetime is flat 
for c < 0. 

It is easy to see that the only spherically symmetric 
MTSs are the spheres given by r = 2M{vq) for a speci- 
fied vq. These will be the apparent horizons on spher- 
ically symmetric Cauchy surfaces which intersect the 
r = 2M{v) surface. Let us denote the r = 2M{v) sur- 
face by H. It is easy to show that H is spacelike and is a 
trapping horizon. The EH lies outside H and is strictly 
separated from H when M > 0; at late times, H asymp- 
totes to the EH 110]. 

Let us now consider non-spherically symmetric 
Cauchy surfaces. There is now an important qualitative 
difference from the Schwarzschild case. There, the ana- 
log of H was null and expansion free; the intersection 
of any spacelike surface with H was then a marginal 
surface, as long as this intersection was, topologically, 
a complete sphere. This is also true more generally 
when the black hole is isolated in an otherwise dynam- 
ical spacetime (if H is an isolated horizon). However, 
in genuinely dynamical situations, H is spacelike as in 
the Vaidya example. In this case, if the intersection of a 
Cauchy surface with H is not one of the spherically sym- 
metric marginal surfaces, then the intersection cannot be a 
marginal surface even if it is a complete 2-sphere. This state- 
ment follows directly from Theorem 4.2 of 1 14]. Thus the 
question naturally arises: are there apparent horizons 
on non-symmetric Cauchy surfaces in the Vaidya space- 
time? One would expect there to be Wald-Iyer Cauchy 
surfaces which come arbitrarily close to the singularity 
but which do not contain any marginal surfaces, but we 
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Figure 2: The t-z section of the Vaidya spacetime where t = v — 
r. The tilted lines are (sections of) the axisymmetric surfaces 
given in eq. Jl} for a range of t values. The marginal surfaces 
on these sections are marked by a The dynamical horizon 
H (the r — 2M{v) surface) are the pair of bold curves, and the 
two dashed straight lines are the v = light cone which is the 
boimdary of the flat portion of the spacetime. The singularity 
is the positive t-axis (z = 0, t > 0). 



shall now see that apparent horizons do exist on a large 
class of non-symmetric Cauchy surfaces. 
We choose the mass function 

, , r fort; < 

^^''> ~ iMov^iv^ + W^) fori;>0 

with the constants Mq = 1 and W ~ 1/10. This is a short 
pulse of radiation that forms a black hole with the final 
mass Mq. This mass function is only at u = 0, but our 
results are unchanged qualitatively for other mass func- 
tions with higher differentiability. For this mass func- 
tion the singular point u = 0, r = is locally naked (see 
e.g. 1 15]), but this is not relevant for our purposes. 

We examine the spacetime with a slicing that is only 
axially symmetric. We use a time coordinate t given by 

t = t — az = V — r (1 + acos6) , (4) 

where t = v — ris the standard Vaidya time, and the con- 
stant Oi determines how much the slice is boosted in the 
z direction. We chose a — 10/11. We have also exam- 
ined other more complicated foliations, but the results 
presented below do not change qualitatively. 

The results are shown in figures |2 and |3] Figure |3 
shows the t-z section of the Cauchy surfaces for a range 
of t values, and the distorted MSs on these sections. This 
clearly shows that the MSs extend outside H and can 
also extend into the flat region. Figure |3] shows the MS 
on the f = -0.3 and t = sHces. The MS at f = is 
a future marginally trapped surface, i.e., ®(n) < 0. At 
t = —0.3, the MS extends into the flat portion, and this 
part of S is planar, with = — 0. On the rest 
of the sphere, 0(„) < as expected; this is therefore a 
weakly marginally trapped surface. Furthermore, the 3- 
dimensional world tube obtained by stacking up all the 
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Figure 3: The location of the distorted marginal surface on 
the x-z section of the Cauchy surfaces at f = —0.3,0.0. The 
solid curve is the intersection with the dynamical horizon H, 
the dotted curve is the distorted MS, and the region inside the 
dashed curve is the c < region, i.e., the flat portion of space- 
time. Note that at t = 0.0, the Cauchy surface does not in- 
tersect the flat portion at all, and the distorted MS coincides 
closely with the DH except near its south pole. In both cases, 
the north pole is inside the DH while the south pole is outside. 



MSs turns out to be spacelike; the ones with 0(„) < 
form a dynamical horizon. 

We also look for surfaces with a small non-vanishing 
expansion ©(^j = ±10"'^. These surfaces can be viewed 
as radial deformations of the MS; the outward deforma- 
tion has 0(£) > 0, and < for the inward deforma- 
tion. At f = 0, the inward deformation is strictly trapped 
and the outward deformation is strictly untrapped. At 
t = —0.3, the inward deformation has > in the 
flat region and 0(„) < elsewhere. The outward defor- 
mation has 0(„) < everywhere and is thus strictly un- 
trapped. We nave not been able to find strictly trapped 
surfaces which extend into the flat portion of space- 
time. Sufficiently far in the future, the MSs asymptote 
to the spherically symmetric DH and also come arbitrar- 
ily close to the EH. Finally, there are restrictions on the 
location of trapped surfaces in the presence of a d5mam- 
ical horizon 1 14]. We have verified that these restrictions 
are satisfied. The existence of such distorted MSs was al- 
ready suggested in 1 14], but with no restrictions on 0(„); 
here we have also shown &{^n) 1^ 0- 



Conclusions: We have numerically studied non- 
symmetric trapped surfaces in simple spherically sym- 
metric spacetimes. We have seen that the Wald-Iyer ex- 
ample illustrates the importance of verifying 0(„) < 
for apparent horizons located numerically. In Vaidya, 
we have found trapped surfaces which extend outside 
the usual r = 2M surface H. This shows that H is not 
the boundary of the trapped region. We have also found 
marginal surfaces that extend into the flat region of the 
spacetime. 

The boundary of the trapped region should be spher- 
ically symmetric, since it is an invariantly defined ge- 
ometric quantity in a spherically symmetric spacetime. 
This lends support to Eardley's conjecture that the event 
horizon is the boundary of the trapped region, since the 
EH is the only natural candidate. However, we have not 
found strictly trapped surfaces that extend into the flat 
region of the spacetime, so that the boundary may be 
inside the EH. 

We conclude with some open questions that need to 
be addressed: (i) Is it possible to push the marginal sur- 
faces arbitrarily close to the event horizon, even in the 
flat region? This would verify that the EH is truly the 
boundary of the trapped region, (ii) For asymptotically 
flat spacetimes, the event horizon is the natural candi- 
date for the boundary of the trapped region. What is 
this boundary for non-asymptotically flat spacetimes, 
e.g. asymptotically de Sitter spacetimes where the event 
horizon is not strictly defined? 
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